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Abstract 

An updated and extended analysis of the quark mass dependence of the nucleon's 
axial vector coupling constant qa is presented in comparison with state-of-the-art 
lattice QCD results. Special emphasis is placed on the role of the A(1232) isobar. It 
is pointed out that standard chiral perturbation theory of the pion-nucleon system at 
order p 4 fails to provide an interpolation between the lattice data and the physical 
point. In constrast, a version of chiral effective field theory with explicit inclusion of 
the A(1232) proves to be successful. Detailed error analysis and convergence tests are 
performed. Integrating out the A(1232) as an explicit degree of freedom introduces 
uncontrolled errors for pion masses > 300 MeV. 



1 Introduction 



The axial-vector coupling constant qa of the nucleon represents a benchmark test of our 
ability to extract hadron properties from the QCD Lagrangian. Its empirical value is ac- 
curately determined from neutron /3-decay [1]: gA — 1-267 ± 0.003. On the other side, 
both lattice QCD calculations [21 El H] and chiral effective field theories [HI EJ [7] are making 
progress in describing the quark mass dependence of this nucleon property. 

The present paper updates and extends our previous study in Ref. [6] about chiral ex- 
trapolations for gAi in the continuum and infinite volume limit. We refer to [8} HI |9] for 
chiral perturbation theory analyses of the effects due to the finite spatial extent of the lat- 
tice simulation volume [10]. In Ref. [6] we compared two different two-flavor, non-relativistic 
chiral effective field theories at leading-one-loop level: Heavy Baryon Chiral Perturbation 
Theory (HB;\PT), with pions and nucleons as active degrees of freedom, and the so-called 
Small Scale Expansion (SSE), which includes the A (1232) explicitly. Treating the A (1232) 
as an explicit degree of freedom turned out to be crucial in order to obtain a consistent 
extrapolation of available lattice data down to the region of small quark masses. 

The important role played by the A (1232) in the physics behind gA does not come as 
a surprise. It has in fact been known for decades that the A-dominance of P-wave pion- 
nucleon scattering, or equivalently, the strong spin-isospin polarizability of the nucleon, 
has its pronounced impact on matrix elements of the axial current in the nucleon ground 
state. The Adler-Weisberger sum rule is a prominent example illustrating this connection. 
Extended versions of chiral effective field theories, such as SSE, have been designed to 
incorporate these well established features. 

A surprising outcome of lattice QCD results for gA is their weak dependence on the quark 
mass. This suggests a subtle balance between contributions from the different degrees of 
freedom involved. In fact, substantial cancellations between irN and 7rA loops in gA have 
been reported first by Jenkins and Manohar [TT] . In Section [2] we demonstrate that no 
interpolation between physical point and state-of-the-art lattice results is possible, with 
parameters consistent with hadron phenomenology, if we use standard chiral perturbation 
theory with only pions and nucleons up to next-to-leading-one-loop order. This observation 
holds both for the non-relativistic and infrared regularized, manifestly covariant scheme 
[12]. Section [3] recalls the Adler-Weisberger sum rule and outlines the importance of the 
A (1232) for in a simple schematic model. In Section I4TT1 we present a detailed statistical 
analysis showing that the leading-one-loop SSE expression with explicit A (1232) worked 
out in Ref. [6] does lead to successful interpolations with recent lattice results even for 
relatively large values of m^. Remarkably, the outcome of our fits turns out to be consistent 
with available phenomenological information about the low-energy couplings involved. In 
Section 14.21 we explore the mapping between expressions worked out with and without 
explicit A (1232) by examining the separate contributions of different powers of m n in the 
expanded SSE result. Conclusions are drawn in Section [5j 

Part of this study has been anticipated in several conference proceedings, see Refs. [T3~l 

mug. 
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Figure 1: Diagrams contributing to the quark mass dependence of g& up to order p 3 in 
Baryon Chiral Perturbation Theory. The wiggly line denotes an external isovector axial- 
vector field, interacting with a nucleon (solid line). The first graph to the left encodes the 
relevant counterterms. 



Figure 2: Graph contributing to nucleon field renormalization at order p 3 . 



2 g A without explicit A (1232) up to Q(p A ) 

The axial-vector coupling constant gA of the nucleon is defined as the limit of the nucleon 
axial form factor at vanishing momentum transfer Q 2 . It is extracted from the forward 
(on-shell) nucleon matrix element of the isovector axial- vector quark current 97^75 {t 1 /2) q 
where q denotes the (u, d) isospin doublet quark field and r % are Pauli matrices. 

In this section we point out the incapability of standard chiral perturbation theory, with 
only pion and nucleon as explicit degrees of freedom, to deal with the quark mass dependence 
of qa- The formalism has been worked out up to next-to-leading one-loop order, 0(p A ), both 
in HBxPT [IB] and in the manifestly covariant framework employing infrared regularization 
[TT] . For the sake of completeness we summarize in Appendix [B] technical details of the 
calculations with infrared regularization, from which the HB^PT expressions can also be 
extracted. 

The diagrams relevant at leading-one-loop level, 0(p 3 ), are drawn in Figs. [1] and [3 In 
those loop graphs all vertices are extracted from the leading nN Lagrangian 

4i = y*{~t/' 9 ^ + 2c v) ^ * + • • • (1) 

where \l/ is the nucleon field; g° A and denote, respectively, the nucleon axial coupling and 
the pion decay constant in the SU{2) chiral limit, i.e. for vanishing u- and ci-quark masses. 
Here a M = a^r*/2 represents an external isovector axial field. 

The next-to-leading-one-loop calculation involves vertices from the second order pion- 
nucleon Lagrangian: 

4* = 4 c ^ m l ^ + f -j* YVK. + (2) 
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Figure 3: Diagrams contributing to qa at the next-to-leading one-loop level. The triangle 
denotes a vertex appearing in C^ N . The wiggly line represents an external axial vector field. 




Figure 4: Graph contributing to the nucleon Z-factor at order p 4 . The square denotes a 
vertex from C^ N involving ci, c 2 and c 3 . 



The trace is taken over the isospin indices. In writing down the first term of Eq. (T2]) we have 
already used the connection between pion- and u-, d-quark masses given by the Gell-Mann- 
Oakes-Renner relation, neglecting isospin violating effects. 

The 0(p A ) result for gA{ m Tr) is obtained by evaluating the diagrams in Fig. [3] and the 
contributions from wave function (Fig.Hj) and mass renormalization, see Appendix [B] The 

(2) 

former ones contain the pion- nucleon- axial vertex from with the two low-energy con- 
stants C3 and C4, which are known to primarily encode the influence of the A (1232) resonance 
on low-energy pion-nucleon dynamics (see for example the discussion in Ref. |18]). 

The expansion of the next-to-leading one- loop expression, Eq. (145 p . around the chiral 
limit gives 



where 
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Here C(A) and F(X) are effective couplings representing unresolved short- distance dynamics 
and compensating logarithmic scale (A-) dependence. The sum of the terms up to and 



including m 3 coincides with the 0{p A ) expression in dimensionally regularized HB^PT [16J. 
2.1 Numerical analysis 

We analyse leading and next-to-leading one-loop expressions in comparison with an updated 
set of lattice QCD results for the quark mass dependence of g^, provided by the RBCK [2], 
LHPC [3] and QCDSF [1] collaborations. 

Fig. [5] summarizes these sets of data, although they refer to simulations with different 
actions, number of flavors, lattice volumes, spacings and procedures to translate the results 
into physical units. The reasons for the evident discrepancies between different lattice groups 
are not yet fully understood, see Ref. [I]. Among the data in Ref. j2] we have selected those 
produced on the larger lattice, with spatial size L = 2.4 fm. The lowest-m^ LHPC point 
corresponds to L = 3.5 fm, see [3]. Both the LHPC and QCDSF collaborations performed 
full-QCD simulations whereas the RBCK data are quenched. Preliminary, unquenched 
RBCK simulations [19] still have too large error bars to make definite statements in the 
region of interest. 

As already mentioned, we translate the quark mass m^-dependence of into a pion mass 
dependence according to the Gell-Mann-Oakes-Renner relation, the leading order linear 
relation between m 2 and m q in Chiral Perturbation Theory. Recent accurate lattice QCD 
results clearly display this behavior for a wide range of quark masses [201 EH E2J [23], even 
for > 0.6 GeV, for reasons not yet understood in detail. 

In a first step the quark mass dependence of qa in chiral perturbation theory is analysed 
without fitting to the available lattice data. We have produced Monte Carlo bands both from 
Eq. (134|) and Eq. (145 p and from the Heavy Baryon formulae at order p 3 and p A . In doing 
this we eliminate g\ by imposing the physical constraint g J 4(mf ys ) = 1.267. The remain- 
ing parameters are randomly chosen within phenomenologically acceptable ranges. The 
dimension- two low-energy constants Ci, C2, C3 and C4 are constrained, within non-relativistic 
chiral effective field theory for pions and nucleons, by several low-energy irN and NN scat- 
tering studies. Combining results of Refs. [21], [25] and [26] and simply superimposing the 
quoted error bars, we obtain the following ranges, in units of GeV~ : 

ciW-1.0 0.7, c 2 « 3.1... 3.5 , c 3 ^-5.6 3.4, c 4 « 3.4 . . . 3.7 . (5) 

The value of c\ basically drives the pion-nucleon sigma term. For a detailed discussion on 
the value of C3, see Ref. [T5]. The analysis in Refs. [23 [2B] of irN — > tttiN scattering at 
leading-one- loop order in HB%PT limits the combination of couplings C r (A) in Eq. (135]) . 
According to [2H], C r \\ = m^ hys ) = —1.4 ± 1.2 GeV" 2 . For the higher-order couplings 
in the infrared regularized expressions — called F and G in Appendix [B] — we rely on 
"naive" dimensional arguments [22], at a regularization scale A = 1 GeV : _F r (lGeV) = 
(— 1 ... 1) GeV -4 , G r (l GeV) = (—1 ... 1) GeV" 6 . For the nucleon mass in the chiral limit 
we scan the range M = (0.88 . . . 0.89) GeV according to the outcome of the analysis in Refs. 

PEED!. 

None of the Monte Carlo bands for g^, which include uncertainties of the low-energy 
parameters, comes anywhere close to the lattice data in the region of interest. This is 
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demonstrated in Fig.Ofor the case of the 0(p 4 ) Heavy Baryon expression. What numerically 
drives the trend of the quark mass dependence of at order p A is the combination of c 3 
and C4 in the term proportional to m\ in Eq. flU) 0. Recoil corrections to the non-relativistic 
results which one gets in the infrared regularization approach, do not improve the situation. 

We can conclude at this point that chiral perturbation theory of the pion-nucleon system, 
at order p A and without explicitly propagating A (1232), fails in the attempt to provide an 
interpolation of gA between the physical point and lattice QCD results. The price one 
would have to pay for enforcing adjustment of an interpolating curve to the lattice data is 
that the combination 2C4 — C3 of irN low-energy constants must be tuned to a value totally 
inconsistent with the empirical ones. For chiral perturbation theory to be able to make 
contact with present lattice data, compensation of the strong trend must come from 
higher powers of m n . In a recent paper [7] a relatively flat quark mass dependence is shown 
to be possibly achieved at two-loop level in the Heavy Baryon framework. In this approach 
the compensation of the m\ term does not arise from the "double log" piece characteristic of 
the two-loop calculation but from a fine-tuning of unknown effective couplings at fifth order. 
The theoretical uncertainties are reported to be acceptably small only for < 300 MeV. 

A more efficient and physically motivated way to successfully interpolate between lattice 
results and the physical point is to include explicit A (1232) degrees of freedom [6]|§ This 
leads to a whole string of higher powers in m n already at leading-one-loop level and the 
outcome of the fits agrees favourably with available information from low-energy hadron 
phenomenology, see Section 14.11 Intermediate A (1232) contributions are well-known to 
play a crucial role in axial current matrix elements because the axial-vector field induces 
strong isovector iV(l/2 + ) — > A(3/2 + ) transitions [32J. The near degeneracy of the A (1232) 
with the nucleon suggests treating both N and A as explicit degrees of freedom. Before 
focusing on the detailed formalism in Section HI it is useful to recall, in the next section, 
some well-known basic physics which outlines the special relevance of the A (1232) in the 
present context of gA- 

3 Adler-Weisberger sum rule and a schematic model 

The Adler-Weisberger (AW) sum rule [33] combines information from low-energy QCD, 
current algebra and dispersion relations to connect low-energy parameters of the irN system 
(f n and gA) with an integral over the difference of the n + p and n~p total cross-sections. 
Using the Goldberger-Treiman relation and the Weinberg- Tomozawa low-energy theorem 
[34] one obtains^ 

g\ = 1 + r % [ W") - °^M\ + ° > ( 6 ) 

Jo W \ M nJ 

1 The large contribution associated with the m 3 term was already pointed out in Ref. [16j . 
2 An extrapolation for gA(rn n ) using a chiral quark model with A(1232) contributions has been produced 
in Ref. [31]. 

3 See, for example, the derivation in |32j . 
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Figure 5: Data base of lattice QCD results for the pion mass dependence of qa provided by 
the RBCK 0, QCDSF [I] and LHPC [3] collaborations. Also shown is the 0(p A ) result for 
Pa^tt) i n Heavy Baryon xPT, with the physical point included as a constraint. The band 
reflects the uncertainty on the input values of the low-energy constants involved. 



where the integral is taken over the pion momentum, q — \q\ = \Juj 2 — m^, and cu denotes 
the pion energy in the nucleon rest frame. The deviation of qa from 1 is tied to pion-nucleon 
dynamics and spontaneous (and explicit) chiral symmetry breaking. The left- and right- 
hand sides of the sum rule turn out to agree at the percent level using = 92.4 MeV and 
an accurate parameterization of the measured ir^p cross-sections |35j . 

In order to demonstrate the prominent role of the A (1232) in the AW sum rule for 
gA it is instructive to perform a simple schematic model calculation as follows. Consider 
the P-wave 7rN forward scattering amplitudes f a {^) in spin-isospin channels a = (21, 2 J). 
Their contributions to the total irN cross section is 

a a (uj) = -q- Im/ Q (a>) (7) 
according to the optical theorem. Next, introduce the i^-matrix K a in each channel by 

1 — i \q | K a {uj) 

K a has poles on the real w-axis located at the physical masses of the corresponding iV and 
A intermediate states. 

Starting from the ttNN vertex Eq. ([T|) and the leading-order 7riVA transition Lagrangian 

4^A = -^^^^ + h. c. , (9) 

J-7T 
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Figure 6: P-wave contributions from the nucleon and A (1232) pole graphs to the total ir + p 
and ir~p cross sections, plotted against the pion momentum in the laboratory frame. In the 
left panel the direct and crossed A (1232)-pole graphs are included, in the right panel not. 



where ty 1 is the Rarita-Schwinger field of the A, the i^-matrix pieces involving direct and 
crossed N and A pole terms are, in the non-relativistic limit [36] : 



K: 



33 



#11 
#13 





~9\ 


12VT/2 


UJ 


q 2 


\ 9a 



+ 



C A 



( A 



^f 2 
#31 = 



A-uj 9A+w 



A c\ 



_ 2uj 9 A 



UJ 



(10) 



where A = Ma — Mn is the delta-nucleon mass splitting deduced from the position of the 
resonance pole in the complex center-of-mass energy plane. 

We use f n = 92.4 MeV on the right-hand sides of these last equations, together with 
A = 271.1 MeV, the delta pole position determined empirically from magnetic dipole and 
electric quadrupole transition amplitudes [37], and ca — 1.5 from the A — > ttN decay width 
(see for example Ref. [IB])- Then the P-wave ttN scattering volumes deduced from Eq. (llOp 
at uj = m n agree very well with experiment, see also Ref. [36]. Eq. (jSJ) and the optical 
theorem lead to the P-wave contributions to the total ir + p and 7r~p cross-sections: 



An 

cr w+p = Im [2/ 33 + / 31 ] 



a 



An 



Im[2/ 3 3 + /3i+4/ 13 + 2/ 11 ] 



(11) 



They are drawn in the left panel of Fig. [61 The curves are in good agreement with the 
empirical ones from near threshold up to q ~ 0.5 GeV [3"5*] I3B] , while the region between 
0.5 and 2 GeV receives contributions from higher resonances. Using the cross-sections in 
Eq. ([6]), we obtain = 1.22. The integral on the right-hand side of the Adler-Weisberger 
sum rule is indeed dominated by the contribution of the A (1232) to the ttN scattering 
amplitude. If the calculation includes only the nucleon Born terms from Eqs. (11 01) . setting 
ca = 0, <J n ± p change as shown in the right panel of Fig. [61 "Switching off" the 7riVA coupling, 
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Figure 7: g& in the Small Scale Expansion at leading-one-loop order. The wiggly line denotes 
an external isovector axial-vector field. All vertices shown here appear in the leading ttN 
and irNA Lagrangians. 
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Figure 8: Graph contributing to wave function renormalization at order e 3 . 



the contribution from the dispersion integral in the Adler-Weisberger sum rule changes sign 
and gA is reduced to 0.99. 

This non-relativistic, tree- level calculation illustrates the role of the A (1232) at the 
physical pion mass. In our previous paper [6] we have also shown that treating the A(1232) 
as an explicit degree of freedom leads to successful chiral extrapolation for pion masses well 
above m5J hys . The lattice results used there are now outdated by more recent and improved 
sets. In the next section we check whether our conclusions in Ref. [6] remain unaltered 
for the most recent input data and perform a detailed statistical analysis, using the same 
methods as described in Ref. [T8] . 



4 g A with explicit A (1232) 

In the framework of the so-called Small Scale Expansion (SSE) [38], the delta-nucleon mass 
splitting A is treated as a small parameter and included in the power counting of "small 
scales" generically denoted by e. The leading one- loop, order e 3 , contribution to g^ is 
represented by the diagrams in Figs|7] and [8l We recall here the result in the continuum 
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and infinite volume case with dimensional regularization [BJ: 
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R= — + A — r-l. 



, 25 



(13) 



is a combination of renormalized third-order couplings, cf. Eq. ( 1351) . and is 
the axial delta-delta coupling; f n and A are the 577(2) chiral limit values of the pion 
decay constant and the delta nucleon mass splitting, respectively. See Ref. [6J for further 
details. The analytic continuation of the previous expression for m n > A is achieved via 
the replacement 




, / A 

-yrril — A 2 arccos 



(14) 



Note that in Eq. ( U2l decoupling of the delta has been implemented up to working order 
[HE]. As a result, g A is the same coupling as in HB^PT calculations without explicit delta. 
This is different in Ref. [9]. 



4.1 Fit results 

We have performed fits based on Eq. ffT2l) and its analytic continuation for m n > A both 
to the two-flavor, quenched RBCK data and the full-QCD, 2 (lighter) +1 (heavier)-flavor 
LHPC results. From the RBCK data we have selected the three points with the lightest 
pion masses. The values of ca and A have been fixed as input. Without any loss of 
generality we have set the regularization scale A equal to 1 GeV. The effective coupling 
C SSE (A = 1 GeV) has been eliminated by imposing that the fit curves pass through the 
physical point. 

Fig. [9] and Table [1] summarize our results for the LHPC data. In drawing the 1-a error 
band (which takes into account correlations between the parameters), we set /„- equal to 
its physical value and ca = 1-5. Furthermore, A = 271.1 MeV from the real part of the 
complex delta pole in the total center-of-mass energy plane |37j . 
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Figure 9: 0(e 3 ) SSE best-fit curves, 68%-statistical error bands and systematic envelopes 
for the LHPC lattice results [3] (left) and RBCK data [2] (right). The physical point is 
included as a constraint. See Table [Q 



The "systematic band" [18] in Fig. [9] quantifies the sensitivity to variations of the input 
parameters. It is given by the envelope of the 1-a bands scanning also the additional input 
values f„ = 86.2 MeV [32], c A = 1.125 [6] and A = 293 MeV (from the 90° nN phase-shift 
in the spin-3/2 isospin-3/2 channel [I]). We point out that the constant g\ appears only in 
combinations multiplied by c\. This is evident from diagram (8) of Fig.[71 In the fit using 
ca = 1.125, the apparent strong deviation of g± from the 577(4) quark model prediction 
9i — $9a/5 ~ 2.2 is therefore of little relevance. The range accessible to c 2 A gi is indeed 
much smaller than that for gi itself, see Table [TJ 

Our estimate of C SSE (1 GeV) is consistent with (limited) information from ttN — > mrN 
scattering. One can indeed link C HB (A) in the Heavy Baryon wN effective field theory and 
C SSE (A) in the framework with explicit delta. Comparing the terms proportional to m 2 . in 
the chiral expansion of Eq. (fT2l) and Eq. (1391 . one gets: 



C HB (A) = C SSE (A) + 



r 2 
( A 



72 TT 2 /^ 



1 23 
5 — Oi 
27 y 



7 



:9a 



25 



9i -9a) ln 



2A 
T 



(15) 



According to our fits we obtain C HB (A = mf ys ) = (—0.45 ± 0.05) GeV -2 , consistent with 
the broad range of values extracted from Ref . [2B] • In agreement with [B] , we stress that the 
"chiral log" in the leading non-analytic quark-mass term is only visible for pion masses well 
below the physical point. 

We have also performed fits to the three LHPC points below 600 MeV in pion mass 
only. The outcome is compatible, within large error bars, with our previous results for 

< 760 MeV in Table ffl 
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Table 1: Fit results for the C(e 3 ) SSE interpolation between the physical value of qa and 
the LHPC points [3]. 





(a) statistical 


error 


(b) systematic envelope 


9a 




1.224 ± 0.004 


fitted 


1.191 . . 


. 1.228 


fitted 


9i 




2.80 ± 0.12 


fitted 


2.6 . . 


. 6.0 


fitted 


C SSE (lGeV) 


(GeV~ 2 ) 


-1.65 ± 0.17 


elim. 


-4.3 . . 


. -1.4 


elim. 


u 


(MeV) 


92.4 


fixed 


86.2 . . 


. 92.4 


scanned 


ca 




1.5 


fixed 


1.12 . 


. 1.5 


scanned 


A 


(MeV) 


271.1 


fixed 


271 . . 


. 293 


scanned 


X 2 /d.o.f. 


0.93 


0.92 . . 


. 1.27 




c a9i 




6.29 ± 0.27 




5.9 . . 


. 7.6 




C HB (mP hys ) 


(GeV~ 2 ) 


-0.45 ± 0.05 




-0.55 . 


. . 0.41 





Table [2] summarizes our study of the RBCK data. The error bars of the output param- 
eters absorb the effects of one heavier flavor, together with the above mentioned systematic 
discrepancies due to different fermion actions, (un) quenching, and translation of lattice 
results into physical units. 

Given the present data situation, extrapolations without the physical point as input 
need to include finite volume dependence in order to get an estimate of qa at the physical 
pion mass with reasonable statistical accuracy. Those studies have been performed by the 
LHPC [3] and the QCDSF collaborations [I], the latter one using our SSE scheme of Ref. [6] 
extended to finite volume. The couplings determined in Ref. jl] agree with our fit results of 
Tabled] Remarkably, both collaborations obtain extrapolated values of qa at the physical 
pion mass that are consistent with phenomenology, although the present level of accuracy 
needs to be improved. 

In the future, reliable chiral extrapolations will certainly be based on combined fits of 
several observables, including finite volume corrections. Such a program will become feasible 
once simulations at smaller pion masses will be available. 



4.2 Expanding in powers of m n 

The effective field theory framework with explicit A (1232) degrees of freedom offers a way 
of studying convergence properties of irN chiral perturbation theory. We perform such a 
test using the 0(e 3 ) non-relativistic SSE expression which is able to describe the quark 
mass dependence of qa over a large range of pion masses. We expand the expression ( [121 
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Table 2: Fit results for the C(e 3 ) SSE interpolation between the physical value of g A and 
the three RBCK points lightest in pion mass with L = 2.4 fm [2]. 





(a) statistical 


error 


(b) systematic envelope 


9a 




1.223 ± 0.007 


fitted 


1.189 .. 


. 1.230 


fitted 


9i 




2.83 ± 0.26 


fitted 


2.5 . . 


. 6.5 


fitted 


C SSE (lGeV) 


(GeV~ 2 ) 


-1.67 ± 0.37 


elim. 


-4.7 . . 


. -1.3 


elim. 


u 


(MeV) 


92.4 


fixed 


86.2 . . 


. 92.4 


scanned 


ca 




1.5 


fixed 


1.12 . 


. 1.5 


scanned 


A 


(MeV) 


271.1 


fixed 


271 . . 


. 293 


scanned 


X 2 /d.o.f. 


0.2 


0.19 . . 


. 0.46 




c a9i 




6.4 ± 0.6 




5.7 . . 


. 8.2 




C HB (mP hys ) 


(GeV~ 2 ) 


-0.43 ± 0.11 




-0.6 . 


. . 0.5 





for m w < A in powers of m w , around the chiral limit. The resulting series is of the form: 

,3 



9A = 9° A <1 



, m n . 2A 
a 2 m — + a 2 In — + p 2 
A A 



A 



2n 



in 



2n 



n>2 



^ n ^2n-2 OA ' / ■> rn /\2n- 



0(e 4 



n>2 



g° A [l + Ami ^ ^ + B(X) m 2 + C ml + Dm% In ^ + E(X) m\ 



(16) 



All terms can be mapped onto the ttN Heavy Baryon expansion, according to the decoupling 
theorem [iQlUT]. The series in m"/A m in the last equation corresponds to "integrating out" 
explicit delta degrees of freedom. The intermediate delta contributions are embedded in a 
string of couplings appearing in the Heavy Baryon Lagrangian. We note that in Eq. (TIB"]) 
the terms starting from m\ come entirely from diagrams 6 - 8 in Fig.0 For m v < A 
a detailed numerical analysis of Eq. (I16p therefore gives an impression of the systematic 
errors introduced by integrating out the leading-order effect of the A (1232) in the HB^PT 
expansion. 

Fig. Hill shows the convergence properties of the series (fT5"|) towards the full 0(e 3 ) result 
in Eq. (1121) . In Fig.[TT]we plot the difference 5g A between the expansion truncated at some 
power of m n and the full expression, for fixed values of the pion mass. In those plots we set 
g A , C ' and g± equal to their central values in the fit to the LHPC data, see Table [H 

For small values of m n the expansion exhibits fast convergence. As soon as m n becomes 
larger than A we are outside the radius of convergence of the series in Eq. (1161) . However, up 
to about 300 MeV in pion mass, the expansion truncated for example at m\. still yields a good 
approximation. This is consistent with Ref. |7] where a fifth-order counterterm proportional 
to m\ turns out to be necessary and sufficient in the Heavy Baryon framework to get in 
contact with the smallest-m T lattice data. 
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m n 2 [GeV 2 ] 

Figure 10: gA in SSE at order e 3 : expansion in powers of m n . Parameters have been set 
equal to the central values in the fit to the LHPC data. The labels n of the curves give the 
order m™ at which the series is truncated. The star denotes the physical point. The LHPC 
data point with smallest m n is also displayed. 



For m n > A the expansion (1161) behaves like an asymptotic series [12] • For a given value 
of m n there is an optimal truncation order giving minimal deviation from the full result. 
Beyond this order any truncation worsens the result. We find also that the optimal order 
shifts to lower powers when increasing the pion mass. At some value of the deviation, 
even at the best truncation order, exceeds the required level of precision. 

Our analysis suggests that for m n > 300 MeV, once the pion mass exceeds the AN mass 
difference, the systematic errors in the Heavy Baryon expansion induced by removing the 
A (1232) as an explicit degree of freedom tend to get out of control. The series in Eq. (jlfip is 
of course only part of the full Heavy Baryon result. It is thus in principle possible that the 
observed deviation 5gA can be compensated by other higher order terms in the diagramatic 
expansion. 

At this point it is worth noting that a reliable expression for the quark mass dependence 
of gA is crucial in order to consistently describe the quark mass dependence of the nucleon 
mass in the Heavy Baryon scheme at order p 5 [131 HI] . Our numerical analysis suggests that 
this cannot be achieved for > 300 MeV in a chiral effective field theory scheme restricted 
to pion and nucleon degrees of freedom only. 

We have explored the sensitivity of our convergence study to g±, which controls the 
relative strength of competing structures in the 0(e 3 ) SSE formula ffl2|) . We observed that 
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for the rather "large" value gi ~ 6 an m^-approximation is quite close to the SSE result for 
a large range of pion masses. We interpret this as a warning that an effective theory can 
spuriously appear to be well behaved at a specific truncation order for a particular artificial 
choice of parameters. This is in agreement with Ref. [7j where the two-loop Heavy Baryon 
expression is found to describe the trend shown by lattice data up to m n m 600 MeV just 
for a particular choice of effective low-energy couplings. 

Our analysis is based on the validity of the SSE perturbative framework. Our fits support 
the conclusion that in SSE important quark-mass dependent effects for g^, involving the 
A(1232), are moved to low perturbative orders in the diagramatic expansion. However, 
before any firm conclusion can be drawn it is mandatory to study the effects of higher-order 
corrections and check the stability of our results. It is needless to say that for a numerical 
study at 0(e A ), we must either wait for an improvement of the statistics for the lattice 
data or perform a combined analysis of g& and other nucleon observables characterized by 
a common subset of low-energy constants. 



5 Conclusions 

Summarizing the main results of this (continuum and infinite volume) analysis of gA, we 
conclude as follows: 

• Heavy Baryon Chiral Perturbation Theory restricted to pion and nucleon degrees of 
freedom only, when applied up to and including next-to-next-to-leading order, fails in 
attempts to produce meaningful interpolations between g& at the physical point and 
present lattice data. 

• The explicit inclusion of the A (1232) at leading one-loop level, as implemented for 
example in the Small Scale Expansion scheme, is crucial in order to get a satisfactory 
description of the quark mass dependence of g& from the chiral limit across the physical 
point up to the lattice data. This does not come as a surprise in view of the well-known 
A-dominance of irN interactions and the Adler-Weisberger sum rule. 

• An expansion of the leading one-loop expression with explicit A (1232) in powers of 
m n exhibits a stable convergence pattern only for pion masses below 300 MeV, the 
characteristic scale of the delta-nucleon mass difference. For larger values of m n the 
expansion behaves like an asymptotic series. Not unexpectedly, our analysis suggests 
that the systematic errors induced in the Heavy Baryon expansion by integrating out 
the A (1232), get out of control for m T > 300 MeV. 

• The addition of one heavier flavor in lattice simulations yields interpolation results 
that are statistically compatible with the two-flavor case. 
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Figure 11: qa in SSE: convergence plot for the expansion of the (9(e 3 ) result in powers of 
m n . The dotted line refers to the full result with parameters set equal to the central values 
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Appendix 



A tf- dimensional one-loop integrals 

We encounter pion loop integrals in d = 4 dimensions of the following basic type: 

1 f d d k 1 



AJm 2 S 



i J (2ii) d m\ — k 2 — ie 
Any ultraviolet divergence as d tends to 4 is subsumed in 

A d ~ 4 f 1 1 



m d ~ 2 (47r)- d / 2 r 1- 



L(A) 



In (4tt) + T'(l) + 1 



16?r 2 [d-A 2 
where A is the dimensional regularization scale. Consequently, 



L(X) 



1 , m„ 
m 



16tt 2 A 



For loop integrals involving a nucleon propagator, consider the infrared singular 
denoted by a subscript / attached to the integral: 



I Nip ,m„ 



d d k 



1 



% J! (27r) d (ml -k 2 - ie) [M 2 - (p - k) 2 - ie) 



p 2 - M 2 + ml 

p2 



L(X) + I N (p 2 J mi) 



I N (p ,m w 



1 ay/1 - n 2 
' 8tt 2 1 + 2aQ + a 2 
1 a(a + Q) 
' 16vr 2 1 + 2aVt + a 2 



arccos 



Q + a 



a- 



with a 



p 2 - m 2 - M 2 
2m n M 



Furthermore, 



p»I { i\p\ml 



d d k 



k'' 



i J j (2vr) d (m 2 - k 2 - ie)[M 2 - (p - k) 2 - ie] 
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Using 

P -k = \{p 2 - Ml + ml) + \{k 2 -ml)- 1 - [(p - kf - M 2 

one finds 



4 1} ml) = ± [(p 2 - Ml + ml) I N (p 2 , ml) + A.(m 2 

We use the following notations: 

I N ee I N (p 2 = Ml, ml), I« = l£V = Ml, ml 

Ia(p 2 ) = In(p 2 , M -h. M°), J A = J A (p 2 = M 2 , m 2 ) 
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B Details on the 0(p 3 ) and 0(p 4 ) calculations 



At order p 3 , nucleon field renormalization contributes in the way shown in Fig.[2j At this 
level of accuracy, the nucleon self-energy S(p) is approximated as 



E(jf = Mat) « E# = M ) , 



(29) 



where Mq is the nucleon mass in the SU(2) chiral limit. This implies for the nucleon 
wave-function renormalization factor: 



dp 1 



--M 



(30) 



Here S^ 3 ^ is the nucleon self-energy at order p 3 . Using infrared regularization [12 



S (3) = |^ (Mo+ ^ m 2 Jjv(p 2 )m 2 ) + (Mo _^^ 7 (l) (p 2 )m 



2- 



(31) 



in terms of the basic integrals in Appendix |A] Hence one finds: 



1 - 



32tt 2 ^M V4 - m2 /M 2 



(2m 3 — QM^m n ) arccos 



m, 



2M n 



/ 777 2 / 770 

+ MjA - ^ (m 2 (1 + 48vr 2 L(A)) - 32n 2 L(X)ml + (3M 2 - 2m 2 ) In ^ 



8fi 20 m 2 +32F 2 m^ 



(32) 



where L(\) is defined in Eq. ( |18p . Here B 2 o is a third-order coupling [6], while _F 2 enters in 
a fifth-order counterterm required to absorb the divergence proportional to m\. 

Projecting out the contributions to qa from the leading-one-loop amplitudes in Fig.fU 
the following expression results, given in terms of the integrals in Appendix [A] 



2m 2 M (d- 3)1 



N 







+ (2(d - 3)M 2 m 2 + ml) g^Mp 



p>=M% 



+ {d- 3)M A n 



32 F 1 ml . (33) 



B 9 is the third-order coupling already introduced in Ref. [6] and Fi takes care of an ultraviolet 
divergence. 

Combining Eq. ( |32l) and Eq. ( 133]) , the result for the pion mass dependence of qa at order 
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p 3 with infrared regularization is given by 



9a 



167rV. 2 M V4-m2/M 2 



9> 3 M9 a + l)M% - (Sgf + 2)ml) arccos ' 



2M n 




Ml 



Mlmlgf + « - 16/ 2 M 2 7T 2 )<& 



+ ((4 9 r + 2)g\mlMl - (3g / + In ^ - 64C-(A)/ 2 M 2 m 2 

+ 32F r (A)mt . 



(34) 



For the couplings involved in third order counterterms we use the notation already employed 
in our previous paper [6]: 

= Bl(\) - 2g A B r w (\) . (35) 

Furthermore, F r (\) = F{(\) + F^X). Those effective couplings are renormalized and 
regularization-scale-A-dependent. They encode short- distance dynamics effects and scale in 
just such a way that the right-hand side of Eq. (134]) is scale independent: 



C r (X) = B 9 - 2g° A B 



L(X) 
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32/ 2 M 2 



-n° + r7° 3 

p \2 

02^ 



9W + Wa1 



(36) 
(37) 



The factor 32 in the fifth-order counterterm in Eq. (1341) emphasizes that F is the effective 
coupling to which one should apply "naive" dimensional arguments, cf. Ref. [2H]. Indeed 
the effective ttN Lagrangian at order p 5 can contribute via 



Cf N = 32fm^a^ 1 ^ + 



where we expect / = 0{\/K 4 



Expanding Eq. (l34|) around m n = leads to 

9°A + 2^ 3 



9 a = 9 A 



+ 



4C r (A) 



03 
9 A 

167T 2 / 2 



2 

T 

3 



ln^ 

Stt 2 / 2 A 



{X) + WftMl + 16vr 2 / 2 M 2 A 



8vr/ 2 M ^ 
1 0(m») . 



711 



(38) 



(39) 



The sum of the first two terms in this expansion coincides with the leading-one-loop 
expression for g A in HB^PT, as expected in infrared regularization. The calculation gives 
a full tower of "recoil corrections" in the form of increasing powers of 1/Mq. We include 
the counterterm 32 F in the third-order calculation in order to achieve renormalization 
without neglecting recoil corrections. Since contact terms up to and including 0(p 3 ) 
cannot absorb higher-order divergences at their /3-functions cannot compensate for 
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N N 
♦ - 



Figure 12: ci-insertion in a nucleon line: the diamond corresponds to the vertex i^c\m\ 
from c!$f. 



scale dependence which is suppressed by two powers of 1/M . By introducing the term 
32F r (A) m 4 , we remove this unphysical scale dependence. 

At order p A nucleon field renormalization (Fig.Hj) and mass renormalization have been 
evaluated consistently within the accuracy in the perturbative, diagrammatic expansion at 
which we are working. In Fig. [12] we draw as a diamond the second order insertion in the 
nucleon line proportional to the low-energy constant c±, see Eq. (EJ. At order p A we have 
to compute all graphs resulting from the insertion of at most one ci-vertex in at most one 
nucleon line in the diagrams (2), (3) and (4) of Fig. [1] and in the diagram of Fig.[2j Since 
for small 



j> - (M - 4 dml) ft — M ft - M 



(iAciml) 



ft- M 



(40) 



we can summarize the ci-insertions by a simple shift of the pole of the nucleon propagator 
from the "bare" nucleon mass to its renormalized value at second chiral order: 



M, 



N 



M -4 Cl ml + O(p 3 ) . 



(41) 



Therefore nucleon field renormalization is taken into account by 



9°A 



1 + l E(3) 



+ If 



(42) 



Here S^ 4 ^(^) is the contribution to the nucleon self-energy from the fourth-order tadpole in 
Fig.H cf. [12], 

^2 



E( 4) = 3^ ( 
J - 



Mid 



c 2 - c 3 



(43) 



What is relevant at order p 4 is just the linear term in the expansion of Eq. (142]) in powers 
of C\. 

Each of the fourth-order diagrams in Fig. [3] contributes to as follows: 



9° A 



(d-l)f2M n 



c 3 [A n + {mi - 4M 2 ) I N ) + c 4 A, + 4(d - 2)M 2 I N + m 2 I 



in 



(44) 
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The pion mass dependence of g A at order p 4 is finally given by 



g A =& + * C r (\) ml + 32 F r (\) ml + 128 G r (X) m 6 n 



32 Cl M$m 4 n - \2c 1 m% 
(4M 2 - ml) Ml 



+ In 



"M 9 A m l 



X J 16vT 2 / 2 L 



167T 2 / 2 

2(1 + 2gf) - (3c 2 + 4c 3 - 4c 4 )^ + (2 + 3^ 2 )^| 



+ ^(24 Cl + C3 + C4 + 36 Cl ^ 2 )^ 



+ arccos — 



m Tl 



n s 



2M ; 16vr 2 / 2 (4M 2 - m 2 ) 3 / 2 



128 



(c 3 - 2c 4 )M 3 



+ 32(1 + <?° )(M 2 + 4 Cl M m 2 ) + 32(c 3 - c 4 )M m 2 - 4(4 + 5g° A )m\ 

4 4 6 

9 777 o 777 o ?77, 

- 8(12 Cl + c 3 + 18c i9a )^ + (2 + 3g° A )^ + 2(24 Cl + c 3 + c 4 + 36c x <& )^ 



+ 



4 



16tt 2 / 2 M 3 



— (9c 2 + 32c 3 - 16c 4 )M 2 - M - ^(c 3 + c 4 )m 2 



(45) 



G r '(A) is a seventh-order effective coupling appearing in the counterterm required to absorb 
the divergence at m\. This compensates the unphysical scale dependence at this power in 
m n . The factor 128 is motivated by the fact that the effective irN Lagrangian at order p 
can contribute via 



^ T 



C% = 128^m^-a^ 75 * + ... 
"Naive" dimensional arguments suggest g = O (l/A^). At order p 4 one finds 



(46) 



F r {\) = F 
G r {\) = G 



32/ 2 M 2 



(3g° A + 2g° A - 3c 2 g A M - 4c 3 g° A M + Ac 4 g A M 



(24c l9A 3 + 16c l9A+ 2 -c 3 gU 



(47) 
(48) 
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